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The Becchi-Rouet-Stora-Tyutin (BRST) supersymmetry is a powerful tool for the calculation 
of the complexity of metastable states in glassy systems, and it is particularly useful to uncover 
the relationships between complexity and standard thermodynamics. In this work we compute the 
Thouless- Anderson-Palmer complexity of the Sherrington-Kirkpatrick model at the quenched level, 
OO ' by using the BRST supersymmetry. We show that the complexity calculated at K steps of replica 

symmetry breaking is strictly related to the static free energy at A" 4- 1 steps of replica symmetry 
breaking. The supersymmetry therefore provides a prescription to obtain the complexity of the 
TAP states from the standard static free energy, even in models which are solved by more than one 
, step of replica symmetry breaking. This recipe states that the complexity is given by the Legendre 

. transform of the static free energy, where the Legendre parameter is the largest replica symmetry 

breaking point of the overlap matrix. 
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■ I. INTRODUCTION 

fi 

' Glassy systems in their low temperature phase always display a very complicated structure of metastable states. 
(— I In fact, this feature can be adopted more generally as the very definition of a complex system, that is a system with a 
Q highly nontrivial structure of its energy, or free energy landscape. Of course, this definition is not limited to systems 
O , with an energy or free energy function, but can be given in general for any system described by a global state function, 
whatever this function is: a fitness function in biology, a cost function in optimization problems, and so on. A common 
' feature to all complex systems, is that the number of metastable states is exponentially large in the size of the system, 
^ in such a way that it is possible to define an entropy density E of the metastable states, normally called complexity. 
lO The calculation of the complexity is in general a crucial task for understanding both the static and dynamical features 
of complex systems. 

In the context of spin-glasses, the problem of calculating the complexity is quite an old one 
\ particular, starting with the classic paper by Bray and Moore a large number of investigations have focused on the 

■ metastable states of the Thouless-Anderson-Palmer (TAP) free energy in the Sherrington-Kirkpatrick (SK) model 
I 0. Despite all these studies, there are still some open questions regarding the TAP complexity in the SK model, 

"j—; ■ mostly related to the consistency of the TAP approach with the standard static calculation of the thermodynamic 
' free energy. A first question, addressed by De Dominicis and Young in 6] , is whether the partition function computed 
^ , within the TAP approach coincides with the standard thermodynamic one. The result of such a study was that the 
I ' two partition functions are the same only if one imposes some identities, whose origin was unclear at the time. On 
■ O the other hand, in P| and 't'I the formal relationships between the calculation of the TAP complexity and that of the 
static free energy were investigated. These studies clearly pointed out that these two quantities were closely related, 
but were unable to establish an exact formal connection between them. In particular, the equations involved in the 
calculation of the complexity are practically intractable in the SK model, where a full replica symmetry breaking 
. ^ solution must be adopted. 

, In a recent work, the problem of the SK complexity has been reconsidered by using a supersymmetric approach |lOj | . 

■ More precisely, it was noted in that the effective action used to compute the number of TAP states with fixed free 
energy density / is invariant under a generalization of the BRST |l2l Il3l Il4| supersymmetry. This invariance can be 
used to generate a set of Ward identities which significantly reduces the multiplicity of solutions of the saddle point 
equations, thus simplifying considerably the calculation of the complexity S. The first important result of ,10j was 
to show that the identities used by De Dominicis and Young to prove the consistency of TAP and standard partition 
functions, were in fact the BRST Ward identities. Moreover, this supersymmetric approach has been used in 

to compute the annealed complexity of the SK model, which turned out to be exactly connected to the quenched 
static free energy at the one step of replica symmetry breaking (1-RSB). This result suggested the existence of a deep 
connection between TAP complexity and static free energy, irrespective of the degree of approximation we use to 
compute them. 



o 
o 



2 



What we show in the present work, is that indeed such a connection is present. The quenched supersymmetric 
calculation of the TAP complexity in the SK model turns out to be completely equivalent to the calculation of the 
standard free energy, with some peculiar connections between the replica symmetry breaking structure of the overlap 
matrices in the two approaches. This proves at the deepest level the equivalence of the TAP and static approaches in 
the SK model, suggesting that such an equivalence may be valid in any glassy system, at least at the mean-field level. 

There is a further motivation to reconsider the connections between complexity and static free energy. In 
a novel method to compute the complexity was introduced, which is independent from the TAP approach, and thus 
is more viable to be used in system where a mean-field TAP free energy cannot be defined. Of course, a key issue 
is whether this alternative definition of the complexity is in general equal to the TAP complexity in those systems 
where a TAP free energy exists. According to this alternative method, in the formulation of [31 , the complexity is 
given by the Legendre transform of the static free energy of r systems forced to be in the same state. This unusual 
free energy has been put in connection in the past with some "replicated" versions of the standard static free energy 
[TtL ITsjl and this Legendre complexity has been thus computed and compared with the TAP complexity. The results 
have been immediately clear in 1-RSB systems, where the two complexities clearly coincide However, this issue 
was up to now far less clear in fuU-RSB systems, as the SK model: due to the nontrivial form of static the overlap 
matrix, it is not obvious what is the correct form of the replicated free energy that has to be used to compute the 
complexity by means of the Legendre transform method. 

The study we perform in this work shows that supersymmetric TAP complexity and Legendre complexity coincide 
in the SK model at any level of replica symmetry breaking. More specifically, by using the BRST supersymmetry 
we show analytically that in order to obtain the complexity at the k RSB level, we have to perform the Legendre 
transform of a standard static free energy calculated at the k + 1 RSB level. The Legendre parameter is the k + 1-th 
symmetry breaking point, that is the largest breaking point of the static overlap matrix Qab- This allows a calculation 
of the complexity starting from the fuU-RSB form of the static free energy of the SK model. Work in this direction 
is in progress 2fy|. 

In section 2 we briefly review the two methods for computing the complexity, and show how, in the TAP context, 
they lead to the same result. Both methods formally require the calculation of the free energy of r systems forced to 
be in the same TAP state, and this is thus the quantity we calculate in section 3 by using the BRST supersymmetry. 
In section 4 we finally show how the quantity we have obtained is related to the standard free energy of the SK model, 
thus giving a general prescription to compute the complexity starting from the static free energy within the Legendre 
transform method. In section 4.B can be found the main original result of our work, namely the formal connection 
between complexity at k RSB level and static free energy at fc -I- 1 RSB level. Conclusions are discussed in section 5. 
A shorter account of our results can be found in 1271. 



II. DIFFERENT METHODS, SAME COMPLEXITY 

In this section we will briefly review the two different methods to compute the complexity, and discuss their mutual 
connections. The first method is only defined when a mean-field free energy, function of the local magnetizations 
rrii, is defined. This quantity is known in the spin-glass context as TAP free energy Ftap, and its local minima {mf}, 
labeled by a = 1, . . . ,Af are identified with the metastable states of the system. The complexity /) of the TAP 
states with free energy density /, at inverse temperature /3, is defined as, 

/) = ^ log ^ 6m f - f3FTApM] = ]^ log / du e^-f'f , (1) 

£1=1 Q = l 

where a indicates a given metastable TAP state. If we compare this equation with the equivalent one in the past 
investigations of 0,13) 0,0,01, we can see that we have made the change of notation u —u. The technical reason 
for this will be clear later. We hope that this change of notation will not make difficult the comparison with former 
studies. If we define the thermodynamic potential as, 

M 

exp(-/3iVu$) ^ ^e-^"^^^^^'"^) , (2) 

a=l 

we can use the steepest descent method to obtain, 

nPJ)=Puf -Pu^{(i,r) , (3) 
where the parameter u — u{f3, /) is fixed by the equation, 

mu)+u'-^^f. (4) 
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In other words, the TAP complexity can be obtained as the Legendre transform of the effective thermodynamic 
potential <i>(/3,u) with respect to the parameter u, which therefore is the Legendre-conjugate variable of the free 
energy density /. 

A different approach to the calculation of the complexity, which does not a priori rely on the existence of a TAP 
free energy, is the one introduced in 15j. The total equilibrium free energy of a super-system composed by r real 
replicas forced to stay in the same state is given by, 

r ^{f3, r) ^ log 1 df e^'^^M ^-PNrf ^ ^xt^ [rf - TS(/3, /)] . (5) 

Note that ^'(/3, r) is the free energy density per replica of such a super-system. The key idea of this approach is that, 
due to the constraint to stay in the same state, there is no degeneracy of the complexity term in the previous formula, 
and thus by tuning r, we can tune the saddle point over /, spanning the entire free energy spectrum of metastable 
states. More precisely, this amounts to say that we can invert relation ((SJ) and find the complexity as the Legendre 
transform with respect to the parameter r of the thermodynamic potential 5'(/3,r), 

S(/?,/)=/3r/-/3rvI/(/3,r) , (6) 

with, 

nP,r)^r'-^ = f. (7) 

Although the formal similarity between equations and © is clear, it may not be as clear what are the physical 
connections between the thermodynamic potentials $(/?, u) and r). In particular, in a standard static approach 
it is not obvious how to compute the constrained free energy 5*. A possibility is to couple in some way the r real 
replicas in the Hamiltonian, and then let the coupling go to zero. As an alternative, we can simply compute the 
normal free energy of a set of r systems by computing the average replicated partition function, Z™, and then impose 
the constraint on the r replicas at the level of the overlap matrix [l7Lll8| . Of course, this needs some careful breaking 
of the replica symmetry, otherwise we simply have, 

lim - log = lim - log[l + rn log Z] = -/3rF ^ -^r-^lr) . (8) 

n^O n n^O n 

On the other hand, in the TAP context it is straightforward to force the r replicas to be in the same state. This can 
simply be done by computing a partition function where we sum over all the r configurations belonging to the same 
TAP state a, and then sum over just one set of TAP states. Such a partition function can be written as, 

Z{r) ^J2Y1 e^^^("^' ■■■Y. e-^"^"-^ . (9) 

a— 1 (Ti^a (Ti-Ga 

Note that of course Z{r) =/= , unless some nontrivial replica symmetry breaking procedure is adopted. The thermo- 
dynamic potential r) is therefore defined in the TAP context by the formula, 

AT 

exp{-(3Nr *) = Z{r) = ^ ^-firFTArim^) ^ (10) 



a=l 



where we have used the standard TAP relation 



^ g-/3//(<T) ^ g-/3FTAj.(m„) ^ (21) 

By comparing (O and 110|l it is now clear that ^'(/3,r) and <I'(/3, w) have exactly the same formal definition, with 
u — r. This identification justifies our choice to use u rather than —u as it was done in such that we 

simply have, 

r) = $(/?, u) . (12) 

Therefore, in the TAP context the two definitions of the complexity coincide, and the key quantity to compute is 
the thermodynamic potential ^'(/3,r) (or, equivalently, $(/3, u)), that is the Legendre transform of the complexity 
with respect to the free energy density /. This quantity is what we calculate in the next section by using the BRST 
supersymmetry. What we will find is that in the SK model ^'(/3, r) is deeply connected to the standard static free 
energy of the system. 
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III. THE SUPERSYMMETRIC QUENCHED CALCULATION 

In this section we will calculate the potential ^'(/?, r) in the TAP context, by using equation H10() . Wc will follow the 
general method introduced by Bray and Moore in 1] , and we will use the BRST supersymmetry, firstly introduced in 
the TAP context in , and discussed for the SK model in JlOj] . We warn once again the reader that for a comparison 
with the previous calculations of P, 0, 0, 0, E| one has to set r — u, since the potential which was normally 
calculated in the past was ^'(/3, —u), rather than \E'(/3, r). 



A. The calculation 

The TAP free energy for the SK model is given by 

FxApim) = -^^Jijin^mj + ^^(j)o{q,mi) , (13) 



with, 



1 _ 

(^0(91 'Ti) = - log(l — TO^) + TO tanh~ (to)— log2 — (14) 

The variables rrii are the local magnetizations, and q is the self-overlap of the TAP states, 

9=^E"^'' (15) 

i 

while the quenched couphngs J are random variables with Gaussian distribution, 

P(Jy) = ^N/27r exp(-iV4/2) . 
The TAP equations and the Hessian of the free energy are respectively, 

/3diFTAp{m) = Jy TOj +0i((7, TOj) = , 

PdidjFTApim) = -PJtj + 4)2iq,m^)Si-j , 

with, 

(j)i{q,m) = 0^ [1 ^ q)m + t&VLii~^ {m) 

= /32(1 - g) + + 0(l/7V) . (16) 

1 — TO^ 

The term of order 1/iV in (/)2((?,to) will be dropped in what follows. From (|10|1 we have that the quenched definition 
of the potential ^'(/3, r) is given by, 



-/3rM'(/3,r) = - logp(/3, r| J) = logp(/3,r|J)" , (17) 
N Nn 



with N ^ 00 and n — > 0, and where, 

M 

p(/3,r|J) ^ ^e-'3'^^-^-("-) = / Y[dm, 5{d,FTAp{m)) \dei{d,djFTAp{ni))\ e-^-F^APir.) (^g) 

a=l i 

We will perform the calculation following closely the lines of ^^id more generally the methods developed in the 
past for this kind of calculation In particular, the modulus of the Hessian determinant will be dropped. This 

approximation is safe only in the energy/temperature region where minima are dominant with respect to unstable 
saddles. We do not know to what extent this is true in the SK model, but we expect it to be true close to the 
ground state (the lower band edge). Since we are interested in the consistency of the calculation of the complexity 
with the statics at the lower band edge, dropping the modulus should be reasonably safe. However, we stress again 
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that this method (as any other method which disregards the modulus) is not under control if saddle points of the 
TAP free energy are exponentially dominant over stable minima: in that case we are weighting each stationary point 
with the undefined sign of its Hessian determinant, with results which are hard to forecast. For a deeper discussion 
of this point see 20]. After introducing the commuting (Bosonic) fields Xi to implement the delta function, and the 
anti-commuting (Fermionic) fields ipi,ijji for the determinant, we find. 



where the action S is given by, 

S{m,x,'4),ip) = '^XidiFTAp{m) + ^tpitljjdidjFTAp{m) - rFTAp{m) 

i ij 

By averaging p(/3,r|J)" over the disorder we obtain the following effective action, 

V 2 



(19) 



(20) 



0^ 



N 



N 



N 



N 



r^a h 



Pi 

2N 



EE 

ab \ i / 

yE E< 

ab \ i / ab \ i 

T 

E [^fMO + i'fi^tMO - rMO] 



ab 



E-X +E E->n -EIE^^'^' 

j J ab \ I / 

\^ n / N \ f ^ 

' -2-E E<^n E-> 



n N 



(21) 



In order to linearize the quadratic terms we follow the standard method to introduce the Lagrange multipliers, 

i 

i 
i 

<5(L„fciV-^<a;,^) = 





e 










dr"^ 












dt''^ 


e~ 




"27ri 






dl''^ 






27r7 







(22) 



and we integrate over qab, Wab, Tab and Lab- In this way the integrals in Xi and (jpi, ipi) become Gaussian and can be 
performed explicitly, giving, 

(23) 



p{P,u\jr ^ J Vn e^^om+N log fU^ dm- 

where n = {qab, A»^ w''\ Wab, i"', Tab, iafc, Z"''} and, 
So(f^) 



ab 



QabLab + W^afc - T^b + y7ah - 2r QabWab 



(24) 



ab 



ab 



h{qaa,m'^) + E'' 



/-I 

'ab 



h{qbb,m'') + ^ 



be c 

w m 



log [Aet{c^2{qaa.m'')5ab + i"")] + -^"'^""^ 



(25) 
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Thanks to the prefactor N in the exponential, the integral in XTi can be performed with the steepest descent method. 
In this way we can write the quenched thermodynamic potential r) as, 



-l3r'^iB,r) = lim - 

n^O n 



So(0) +log / Wdrn'^e 



where, as usual, Q is solution of the saddle point equations 



with, 



mm))) - 



1 



The following saddle point equations are easily solved: 



dWat 
dTab 
dLab 



1 



Wab = -^w'^'^ + rqab 

P 



1 



-t 



ah 



b 



= ^ = 9- 



In order to have expressions as similar as possible to previous investigations 0, ^3 > '^^ define, 

Bab - /3'(l-gaa)(5ah+<"' 
-A„6 - /32(l-(7aa)5a6+U'"'' , 

and by using the explicit forms of 0i, 4)2 we finally obtain. 



2/3' 



c{n,m"-) 



E 

ab 

—r 



ab 

52 



~^ lab + ^ 9a& ^ 9a6 



--log[(2^/32)"detgafc] 



^ 0o(gaa, m"^) + ^"'"rn^'m'' + logdet (t^^ + Bab) 



—y 

2/32 



ab 



tanh"^TO''-^A° 



<lab 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



This quantity has to be extremized with respect to the variational parameters Qab, Bab, ^ab, ^ab- This task is technically 
very hard if we assume a nontrivial form for the various matrices. In particular, if we want to perform a /e-RSB 
calculation of the complexity, at any value of k, the situation is practically hopeless. What we shall see in the next 
section is that the BRST supersymmetry dramatically simplifies the equations, leaving basically just one matrix, qab, 
to be fixed variationally, as in the static calculation. 



B. Using the BRST the supersymmetry 

In nH it was shown that the action in (|2()(l is invariant under a generalization of the Becchi-Rouet-Stora-Tyutin 
supersymmetry |il2..13j: if e is an infinitesimal Grassmann parameter, (|20|) is invariant under the transformation. 



Sirii = eTpi Sxi — eripi Sipi — —exi Sipi = 



(33) 
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The invariance of the action impHes that also the average of any function of the fields m, ip, V', x must be invariant as 
well. Thus, if we set to zero the variation of mi?/5i and Xiipi \l(X\l^. we obtain the two BRST equations, 

+ {m^Xi) = (34) 

r i'lpiipi) + (XiXi) = . (35) 

Once we replicate the action, from equations (|22() and H30|l we have. 





= Tab = 


~Bab/P^ - 


^ (1 - qaa) Sab 




= Rab — 


-^ab/P^ 


- (1 - Qaa) Sab + qr 




= Lab ■ 







(36) 

The first BRST equation therefore becomes, 

Aab^-Bab + (3\abr . (37) 

In order to use the second BRST relation we need the expression for Lab, which can be obtained by the saddle point 
equation d'So/dqab — 0, where Sq is given in 124|) . By doing this, and by using the first BRST equation, we obtain: 

A"" ^\^ab- (38) 

As we see, thanks to the BRST relations the variational parameters Aaj, and \ab can be expressed as functions of 
Bob-, Qab and the parameter r. From the saddle point equations for the variables A'^'' and Bab we get the equations, 

qab = ((m"m^)) (39) 

Bab = -/32(l-<7aa)<Sa6 + /3'((7j-logdet ( -^^ + Bah ))) (40) 

dBab \l-mi J 

which are sufficient to fix Bab and qab- It is possible to prove that the remaining saddle point equations are automat- 
ically satisfied by the BRST expression for A^h and \ab, and therefore do not need to be considered. By using the 
general formula. 



91ogdetM„f, 



{M'^)ab (41) 



dMab 

it is easy to show that equation H4U|I admits the solution Bab = 0. As in the annealed case, this is the solution we 
adopt. Thus, the two BRST relations we are left with are, 

Aab^P\abr (42) 



\ ab 



1 



-B'^r'^qab , 



(43) 



and the only unknown parameter left is the matrix qab- If we use the two relations (I42|l and H43|) in equations H31|) 



and H32|l . and make the change of variable 



quenched thermodynamic potential ^'(/3, r), which is one of our main results: 



ft," ~ tanh (m"'), we obtain a much simpler expression for the 



/3^'(/3,r) = -log2 



An 



r 

ab 



Qaa) 



/n 
Yldh" 



[27r/32 deiqabf'^ 



with, 



n n 

T{h'';qab,r)^^—Y, h^-^ h'' + rj^ log cosh 



(44) 



(45) 



ab 



In the formula above the limit 

-.a^b 



* is understood, and the matrix qab has still to be fixed by the saddle point 
equation qab = ((m'^m*)), with m" — tanh(ft,'^) and where the distribution p(/i°) — exp[J^(/i")]/(27r/3^ det gab) must 
now be used to compute the average ((•)). 

As we have seen in section 1, the complexity S(/9, /) is just the Legendre transform of '^{(3, r). Thus, the quenched 
supersymmetric TAP complexity of the SK model is given by. 



S(/3,/)=/3r/-/3rvE'(/3,r) , 
where the parameter r = r(f3, f) is fixed by the equation, 

a^'(/?,r) 



vE'(/3,r) 



dr 



= f 



(46) 
(47) 
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C. A special case of the BRST supersymmetry: the Bray-Moore action 

Before proceeding, we discuss here a point which was shghtly confusing when the comparison was made between 
the past calculations of the complexity and the most recent ones. In all the classic calculations (both annealed 
and quenched) performed by Bray and Moore 0], and also by Bray, Moore and Young 0, the J-dependent part 
of FTAp{m) in the action was eliminated by using the equations diFTAp{Tn) = enforced by the i5-function. More 
specifically, it was used the equation. 



1 



\jm^m.j = -^^TOj0i(g,mj) , (48) 

ij i 

which is valid in the TAP states. This substitution simplifies considerably the calculation, but the action obtained 
in this way is no longer invariant under For this reason, in our former calculation of [10|. as in the present one, 

we used the full form of Ftap(p^), equation H13|l . and due to this the comparison of our results with those of 0,0,0 
proved somewhat difficult. Moreover, it was not clear what was the equivalent of the BRST complexity, when one 
used the 'trick' H48I) . which seemed to break the BRST invariance from the outset. 

The situation gets much clearer once we realize that the action used by Bray and Moore is actually invariant, under 
a slightly modified version of the BRST supersymmetry. More precisely, the Bray-Moore action, 

SBM{m,x,i),^^) = '^XidiFTApim) + ^i>itl)jdidjFTApim) -rGBuim) , (49) 

i ij 

where, 

GBM(m) = -^^TOi0i(g,TO,) + ^^(l)o{q,mi) , (50) 

i i 

is invariant under the following modified BRST transformations, 

Snii = ei/ji Sxi = ^er^jji S^.^ = -e ^Xi + ^ rm,j^ = , (51) 

to be compared with H33|l . Choosing the same two observables as in the former section, we get the modified BRST 
identities, 

(V'iV'i) + {miXi) - {m^mi) = (52) 
ir {miXi) + {tpiilJi) + (xiXi) = , (53) 

(we recall that in order to make a comparison with former calculations we have to set r = —u). In terms of the 
variational parameters introduced in the previous sections, and once we set B = 0, we obtain, 

Aab^^l3\abr (54) 

Xab = lp^r\ab . (55) 
It is interesting to observe that if we make the change of variable suggested in [2^ , 

Aab Aab+^P^rqab (56) 

3 

Kb + 8 P'^r'^Qab , (57) 

we can reduce relations (I54|l - (|55|l to (|42|I - (I43I), and accordingly reduce the action of Bray-Moore to our action ijSU- 
(|32|l . This is consistent with the results of [23, and shows that the two calculations performed with or without the 
'trick' l|48|l are connected by a simple change of variables, under which the BRST supersymmetry is conserved. 

The important point is that equations l|54 |l - H55|) exactly coincide with the 'ansatz' used by Bray, Moore and Young 
in their calculation of the quenched complexity at the full-RSB level (relations (19) of 7]), which was consistent with 
the previous Bray-Moore quenched calculation of 0. This fact proves that the quenched complexity considered in 
00 was in fact BRST symmetric. On the other hand, as it was shown in [Tol |. the annealed complexity of the total 
number of TAP states considered by Bray and Moore in was not BRST symmetric. For an extensive discussion of 
the comparison between the annealed BRST complexity of ^of and the BRST-breaking one of , see and . 
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IV. CONNECTIONS WITH THE STATIC FREE ENERGY 



Now that we have obtained a general quenched expression for the thermodynamic potential ^'(/3, r), we want to 
investigate what are the connections with the standard thermodynamic potential, that is the free energy of the system 
F{(3). We recall that the thermodynamic potential \E'(/3,r), is the constrained free energy density (per replica) of r 
real replicas forced to stay in the same metastable state, and thus by definition we have F{fi) = ^{f3,r = 1). This 
identity was first proved by De Dominicis and Young in "g], with the assumptions of a particular ansatz, which was 
shown in to be nothing else than the BRST relations. Our goal now is to investigate further the relation between 

and F at a generic value of r 7^ 1. 



A. Complexity vs statics: a preliminary step 



First of all, we note that the annealed calculation of 5'(/3,r) is equivalent to assume for the TAP overlap matrix in 
equation the simple form Qab ~ q5ab- If this is done, it is straightforward to prove (as it was done in that 
the annealed potential ^' is equal to the quenched static free energy, calculate at the one step of replica symmetry 
breaking, with self overlap qi = q, mutual overlap go = and replica symmetry breaking point x — r. In the annealed 
case we obviously have just one value of the overlap and therefore in we could only find a connection with 
FiRSB, once we set go = 0. Now that we have done the quenched calculation it is natural to expect that the ORSB 
potential ^oi?SB is connected to the IRSB static free energy FmsB, with mutual overlap go 7^ 0. Thus, before 
considering the case of a generic number of steps of replica symmetry breaking, we will focus on this simpler case. 
We assume a ORSB form of the TAP overlap matrix qat in H44(l . that is. 



qab = go + {qi - qo) Sab [ORSB] . 

The annealed case is recovered by setting go = 0. By using this ansatz in equation H44|l . we find 

/?*(/?, r) = - log 2 + ^[(r - l)qj + 2gi - 1 - rgg] - ll{q,,, r) 



(58) 



(59) 



where we have defined, 



I{qab, r) 



qo 



2(gi - go) 



log / l[dh- 



V2^/32(gi-go) ' 



and the limit n ^ in the expression for detgaf, has already been taken. In order to proceed we note that, 

[q^\b = ~qo/{qi - qof + ^abl {qi - qo) 

The term X^afc ^"^lab^^ the exponential in J- can be thus rewritten as. 



(60) 



(61) 



qi - qo 



qo 



(qi - qo^ 



(62) 



By using the following Hubbard-Stratonovich identity 



exp 



90 



2f3Hqi-qoy 



dz 



V27rgo 



exp 



2go /3(gi - go) 



(63) 



we find that replicas factorizes in the integral in /(gat, u) and thus we can pass from the n variables ha, to one single 
scalar variable h. After making the further change of variable h — (3{z + y), we finally obtain. 



I{qab,r) = 



dz 



\/27rgo 



2<J0 lop 



dy 



^27r(gi - go) 



2(ii-<!o) cosh'' [I3{z + y)] 



(64) 



If we substitute this form into equation H59() , we find that we have exactly reconstructed the static free energy of the 
SK model at the IRSB level of approximation [2^ W^ . where the replica symmetry breaking point (normally called 
x) is equal to the parameter r. Therefore we have. 



ORSB 



F 



IRSB 



(65) 
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We stress that the overlap matrix Qafs of the static free energy has a different structure from the overlap matrix qat 
of the potential More precisely, for ^oflss "we have the two parameters overlap matrix, 

Qab = qo + {qi - Qo) ^ab ■ (66) 

Moreover, '^qbsb is a function of the variable r. On the other hand, for the static free energy FmsB we have the 
three parameters overlap matrix Qap, 

QaP = go + (gi - go) ^al + (l - gi) ^al3 , (67) 

where ej^^ is a block ultrametric matrix, equal to 1 within a diagonal block of size r and zero elsewhere. In other 
words, the variable r of the potential ^' calculated at the ORSB level becomes the replica symmetry breaking point of 
the static free energy F calculated at the IRSB level. 

This result is in accordance with what found in other systems solved exactly by one step of replica symmetry 
breaking, where the ORSB complexity is alwa ys r elated to the IRSB static free energy by means of a Legendre 
transform with respect of the breaking point x [l^- In the SK model, however, the statics is solved by a full replica 
symmetry breaking ansatz, with an infinite number of breaking points. It was therefore not clear whether the Legendre 
relation between complexity and free energy was preserved in the SK model. More practically, once we consider /cRSB 
solutions, we have more than one breaking parameter x, and thus it is not obvious which one of the k breaking 
points must be used to perform the Legendre transform. What we show in the next section is that the constrained 
thermodynamic potential ^(/?, r) computed at the k RSB level is identical to the static free energy F(/3) calculated 
at the fc + 1 RSB level, with r being equal to the largest breaking point Xmax- This means that the Legendre relation 
between TAP complexity and static free energy is conserved at any level of replica symmetry breaking, and that the 
Legendre parameter is the largest replica symmetry breaking point. 



B. Complexity vs statics: the general case 



We start by recalling the general from of the static quenched free energy in the SK model, before any ansatz on the 
overlap matrix is done |^] , 



a>/3 



a2 



(68) 



where Qap is a ris x Us matrix, with Us — > 0. The subscript in Ug stands for static^ and it is necessary in order to 
distinguish the size of the static overlap matrix Q^p from the size n of the TAP overlap matrix qab- As a first step to 
make the general expression (|44ll for ^'(/3, r) a bit closer to the static free energy (|68|) . we can write. 



COS. 



such that in the expression of ^' we can integrate now over the Gaussian variables ha and obtain. 



'> - T + 1; <' - '> S + 1;^ S + li S - i ">« S 

a a>b a fr*'! 



n2 " r 
ab tJ.v 



T^QabTb 



(69) 



(70) 



We want to prove that ^'(/3, r) = F{(3) if some suitable structures for the overlap matrices are considered. It is 
straightforward to see that this equation is fulfilled if the two following relations are satisfied. 



1 



— log V exp 

nr ^ — ^ 



Y I E E ^alabT^ - E E 
\ ab liv 



a 

r 



— log V 



exp 



fD2 » = 



T — 1 \ ^ 2 r \ ^ 2 1 \ ^ 2 

■2^1aa + -2^<lab = —l^Q^P 

a a>b a>P 



2n 



The first thing to note is that the form of the equations suggests the identity, 

Us = r ■ n . 



(71) 
(72) 

(73) 
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Once this identification is done, we can connect the cr" spin variables {a = l,...,??.^), to the spin variables 
(a = 1, . . . , n; fj, — I, . . . ,r) in the following way, 

{ai,...,(T„J^ {tI,...,t^,t^,...,t^, ,T^,...,0 (74) 

Let us now assume that we are performing the calculation of the potential vE'(/3, r) at k RSB level. In this case the 

(k) 

TAP overlap matrix q^f^ is given by, 

fc+i 

<l^l^-1o + T.('l^-1^-^)'iT\ (75) 

i=l 

with, 

Vk+l = 1 , 4b' = Sat . (76) 

The matrices are the n x n ultrametric block matrices, equal to one on the diagonal blocks of size yt and zero 

elsewhere. The variables jji are thus the replica symmetry breaking points. Unlike the standard static case, in the 
TAP approach the diagonal of the overlap matrix qaa — Qk+i is essential, as it contains the self-overlap of the TAP 
states. For this reason the largest breaking point is trivial, that is yk+i = 1- Therefore, there are k + 1 values of the 
overlap, but only k nontrivial breaking points, ans thus Qat is indeed a fcRSB matrix. The simplest example of this 
matrix was shown in the last section. 

By using relation H74|l it is possible to prove the following key formula, 

n r rn 
ab ti,v aP 

and we recall that rn — rig. It is now possible to prove that, 

n r n r rn 

E E -U'^^i^ - E E '^i'^^ = E QiT -"-^ ' (78) 

ab tJ,v a fj, a#/3 



with, 

QiT - +'^('1^ - 'i^-i) ■ (79) 



Equation (|71|l is therefore verified by the rn x rn ultrametric matrix Q^,^ , which has a standard RSB form with 
fc + 1 levels of replica symmetry breaking. The entries of Q^'^^^ are the same as the TAP overlap matrix g^^^ with 

the only important difference that the elements on the diagonal of Q]^^ are irrelevant, since the sums in 1)68(1 are 
only performed for a ^ (3. On the other hand, it is clear from 1(791) that the k + 1 replica symmetry breaking points 
(cci, . . . , Xk+i) of the matrix Q^^f^^^ are given by, 

xi = ryi 
Xk = ryk 

Xk+i = ryk+i = r . (80) 

In other words, the parameter r represents the largest breaking point of the static overlap matrix 
Finally, from the simple relation, 

n rn 
\^ Jn.Vi} _ 1 \^ (rn.ryi) ^ni \ 

a=^b a^p 

it is straightforward to check that also relation 1(72(1 is verified. 
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Summarizing, what we have proved is that the thermodynamic potential 5'(/3,r) calculated at the k RSB level is 
equal to static free energy F{f]) calculated at the fc + 1 RSB level. In other words, the free energy of r real replicas 
forced to be in the same state, computed at the k RSB level, is equal to the ordinary free energy of one single system, 
computed at the k + 1 RSB level, where the extra k + 1-th symmetry breaking point is equal to r. The symmetry 

breaking points of the static matrix Q^^^^^ are simply the ones of the TAP matrix q^^^ rescaled by the parameter r. 
We can write our result as, 

where the relation between q'^ and Q^^"*"^^ is given above, and, as we have seen, the parameter r is the largest replica 

symmetry breaking point of the static matrix Q^^j,^^^ ■ From equation 1)46(1 . and given the relation between r) and 
F{j3), we finally have the general Legendre equation connecting the supersymmetric complexity of the TAP states to 
the standard static free energy in the SK model, 



with X fixed by the Legendre equation. 



n(iJ)=Pxf-pxF{(i-x) , (83) 



f = F{p-x)+x^-^^. (84) 



This result can be summarized by saying that the quenched complexity of the TAP states is the Legendre transform 
of the static free energy with respect to the largest breaking point r of its overlap matrix. Note that from equation 
(|83|l it is trivial to check that consistency with the statics is obtained at any level of RSB. Indeed, S = for f = F, 
i.e. the lower band edge coincides with the static free energy; moreover, the derivative of S with respect to the free 
energy / at the lower band edge is just Px, where x is the largest static replica symmetry breaking point, as it should 
be. The fact that this consistency can be so transparently read from equation H83I) is one of the main virtues of the 
Legendre-supersymmetric approach. 

If we invert our point of view, and fix the largest breaking parameter of the static overlap matrix Qap to its 
equilibrium value Xmaxi and let r free, we can give a deeper interpretation of these results. First, let us recall that 
the symmetry breaking points of the overlap matrix are related to the overlap probability distribution P{q) |24j . In 
particular, the probability of the TAP self-overlap qaa, is given by, 

Waa = 1 - Umax ■ (85) 

Therefore, if we have, 

r ^ Xyyiax Umax — ^ 1 i (^^) 

r 

the probability of the TAP self-overlap, Waa is different from zero, and thus the number of TAP states is not 
exponentially large in N. This is consistent with the general philosophy of the Legendre transform method which 
states that if the parameter r is larger than the static equilibrium breaking point Xmax , we are stuck at the equilibrium 
ground states, and therefore / = F{(3) and from equation 184|l we obviously find E = 0. Thus, the complexity of the 
static equilibrium states is by definition zero. If, on the other hand, 

r <C Xjyiax Umax — ^ 1 : (^''^) 

r 

then the weight of the TAP self-overlap is zero. In this phase of r the number of TAP states is exponentially large 
in N, their complexity is nonzero and their free energy density / is larger than the equilibrium value F{j3). If we 
parametrize the ultramctric matrices qab and Qafi by means of the two functions q{y) and Q{x) respectively, we can 
write, 

Q{r ■ y) = q{y) . (88) 

This last relation explains quite well the fact that the standard static approach and the TAP approach are simply 
connected by a coarse graining relation, where r is the scale parameter: what we are doing by using the TAP 
approach is to start at a higher level in the coarse graining procedure, since the elementary objects are states, rather 
than configurations as in the static approach. We would like to remark that the many of the physical considerations of 
this last paragraph were already pointed out in where, among the other things, it was for the first time uncovered 
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the Legendre transform relationship between complexity and free energy, and the role of r as an exploring tool of the 
spectrum of states. However, the fact that the modified BRST symmetry (|51|l of that calculation was not recognized, 
prevented to prove equation H83|l . and the consistency between TAP complexity and static results was quantitatively 
tested only close to the lower band edge and to the critical temperature. In our case, from (|83ll we automatically have 
this consistency at any value of / and T. 

A final point to note is that if we set r = 1 the two matrices g^^'' and Q^^f^^^ are actually the same, since the largest 
breaking point of Q^f^^'' is one, and thus this is in fact a /cRSB matrix. This is consistent with the fact that for r = 1 
the thermodynamic potential r) is just the standard static free energy [Sj 



V. CONCLUSIONS 



The main result of this work is given by equations H83|l and (|84|l : the BRST quenched complexity of the TAP states 
in the SK model is just the Legendre transform of the static free energy. The key tool for obtaining this result has 
been the supersymmetry. Our result can be reinterpreted by saying that the degree of difficulty of the computation 
of the complexity becomes equal to the one of the standard free energy once the BRST symmetry is used. This is an 
important point. In principle, from a technical point of view these calculations are quite different. For the complexity 
we need the TAP free energy Ftap, its derivative and its Hessian, while for the static partition function we simply 
need the Hamiltonian H. Even though the relationship between Ftap and H may be easy to uncover Q, it still 
would seem that the calculation of the complexity has a higher degree of difficulty, since it also involves BFtap, and 
d^FxAP- This is exactly the technical problem that the old calculations of the TAP complexity had to face, whenever 
a connection with the statics was investigated. These calculations treated Ftap, QFtap and O^Ftap as independent 
functions, while clearly they are not. In so doing an important physical information was wasted. 

The BRST supersymmetry exactly takes care of the fact that Ftap^QFtap and O^Ftap are not independent, 
and therefore it reduces the redundant difficulty of the calculation, making it practically equivalent to the one of 
the standard free energy. This fact suggests that the BRST supersymmetry should be considered an essential tool 
each time the complexity of a glassy system has to be computed, since it encodes at the deepest level the natural 
connections between the state function (which in structural glasses may be the Hamiltonian) and its topological 
properties, expressed in the distribution of the metastable states. However, we know that in general there may be 
some solutions of the saddle points equations which break the BRST symmetry, still giving a nontrivial complexity. 
At the moment, it is unclear what is the physical meaning of these non-BRST saddle points, nor whether they should 
be preferred to the BRST ones [25112^ . What we believe it can be said with some confidence, is that the consistency 
with the statics, i.e. the equalities between static ground state and lower band edge, and between breaking parameter 
X and is given by the BRST saddle point. Indeed, as already stressed, the present method is only valid in the 

energy phase where stable minima are dominant over unstable saddles, and this is certainly true at, or very close to, 
the lower band edge. Whether the BRST solution of the saddle point equations is the only relevant one also at higher 
energies, depends on the mutual distribution of minima and saddles. Of course, the same uncertainty also holds for 
any other non-BRST solution. 

As we have seen the complexity of the TAP states is given by the Legendre transform of the free energy with 
respect to the largest breaking point of the overlap matrix. It has been recently conjectured in p9l | that when there 
is more than one step of replica symmetry breaking the complexity of the clusters at level i may be given by the 
Legendre transform with respect to the breaking point Xi. This seems a sound generalization of our result, since the 
deepest clusters, associated to the largest breaking point, are indeed the states. Moreover, this conjecture raises the 
interesting issue whether, in some cases, it may be more relevant from a dynamical point of view to calculate the 
complexity of the clusters, rather than the one of the states. It would be interesting to study whether the conjecture 
of can be exactly proved, perhaps using the supersymmetric approach of the present work. 

We thank Alan Bray, Luca Leuzzi, Mike Moore, and Tommaso Rizzo for many stimulating discussions. 
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